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THE TWO-TYPE RICHARDSON MODEL WITH UNBOUNDED 
INITIAL CONFIGURATIONS 
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The two-type Richardson model describes the growth of two com- 
peting infections on Z'^ and the main question is whether both in- 
fection types can simultaneously grow to occupy infinite parts of Z'*. 
For bounded initial configurations, this has been thoroughly studied. 
In this paper, an unbounded initial configuration consisting of points 
X — {xi, . . . ,Xd) in the hyperplane Ti. — {x £ Z'' : xi = 0} is consid- 
ered. It is shown that, starting from a configuration where all points 
in 7i\{0} are type 1 infected and the origin is type 2 infected, there 
is a positive probability for the type 2 infection to grow unboundedly 
if and only if it has a strictly larger intensity than the type 1 infection. 
If, instead, the initial type 1 infection is restricted to the negative xi- 
axis, it is shown that the type 2 infection at the origin can also grow 
unboundedly when the infection types have the same intensity. 

1. Introduction. One of the simplest models for spatial growth and com- 
petition is the Richardson model, introduced in Richardson [13]. The origi- 
nal version describes the growth of a single infectious entity on Z'^, but the 
mechanism can be extended to comprise two entities, making it a model for 
competition on Z"^; see Haggstrom and Pemantle [7, 8]. This paper is con- 
cerned with the two-type version of the model in d > 2 dimensions, started 
from a configuration where one of the entities occupies one single site in an 
infinite "sea" of the other entity. 

The dynamics of the one-type Richardson model are such that an unin- 
fected site becomes infected at a rate proportional to the number of infected 
nearest neighbors and, once infected, it never recovers. This is equivalent 
to first-passage percolation with i.i.d. exponential passage times. The main 
result, dating back to Richardson [13] and Kesten [10], is that, for bounded 
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initial configurations, the infection grows linearly in time in each direction 
and, on the scale the set of infected points at time t converges almost 
surely to a deterministic shape; see Theorem 2.1. 

In the two-type version of the model, a second infection type is introduced 
and the two infections, referred to as type 1 and type 2, respectively, grow 
simultaneously on Z,'^, the dynamics being that an uninfected site becomes 
type i infected at a rate proportional to the number of type i infected near- 
est neighbors and then stays (type i) infected forever. The model has two 
parameters, denoted by Ai and A2, indicating the intensities of the infections. 

For the two-type Richardson model, the main challenge lies in deciding 
if it is possible for both infection types to simultaneously grow to occupy 
infinite parts of Z*^. We will denote by Gi the event that the type i infection 
reaches sites arbitrarily far away from the origin. It has been conjectured 
that in the two- type Richardson model in d > 2 dimensions, started from one 
single site of each infection type, the event Gi D G2 has positive probability 
if and only if Ai = A2; see Haggstrom and Pemantle [7]. The if direction of 
the conjecture was proved for d = 2 by Haggstrom and Pemantle [7] and, 
in the general case, independently by Garet and Marchand [6] and Hoffman 
[9]. 

As for the "only if" direction of the conjecture, the best result to date is 
that, with one of the intensities fixed, mutual unbounded growth has proba- 
bility for all but at most countably many values of the other intensity; see 
Haggstrom and Pemantle [8] . The full result would indeed follow from this if 
it could be proven that the probability of the event Gi H G2 is monotone, in 
the sense that it decreases as the difference between the intensities increases. 
This monotonicity is certainly believed to hold on the Z'^-lattice. However, it 
has turned out to be very hard to prove and, in fact, there are other graphs 
where it actually fails; see Deijfen and Haggstrom [2]. On the other hand, in 
Deijfen and Haggstrom [3], it is shown that, as long as we restrict ourselves 
to bounded initial configurations on Z*^, the particular choice of initial sets 
is not important in deciding whether the event of mutual unbounded growth 
for the two infection types has positive probability or not. 

The purpose of the present paper is to study the two-type Richardson 
model with unbounded initial configurations. The question is as before: can 
both types simultaneously infect infinitely many sites? With infinite initial 
configurations for both types, the answer is (other than in absurd cases) 
obviously "yes." We therefore restrict to cases where type 1 starts with 
infinitely many sites and type 2 with finitely many (for concreteness, only 
one) and the question of coexistence then reduces to that of whether type 
2 can survive. Write {xi, . . . ,Xd) for the coordinates of a point x S Z'^ and 
define TC = {x : xi = 0} and C = {x : xi < and Xi = for all i>2}. Writing 
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for the origin, the initial configurations that we wiU consider are 

ah points in TC\{0} are type 1 infected and 
is type 2 infected; 

ah points in >C\{0} are type 1 infected and 
is type 2 infected. 

Interestingly, it turns out that the set of parameter values that allows for 
coexistence is slightly different for these two configurations. Let P^q^ (resp. 

q"^^) denote the probability measure associated with a two-type process 
started from configuration I {71) (resp. /(£)) and note that, by time scaling, 
we may restrict our attention to the case Ai = 1. Our main result is as 
follows. 

Theorem 1.1. For the Richardson model in d>2 dimensions, we have: 

(a) PjI\{G2) > if and only if X> 1; 

(b) Pc,o(.G2) >0 if and only if \>l. 

In words, when the type 2 infection is strictly stronger than the type 1 
infection, there is a positive probability for the type 2 infection to grow un- 
boundedly in both configurations. Intuitively, the type 2 infection can use 
its higher intensity to rush away along the xi-axis and achieve an unsur- 
passable lead over the type 1 infection. When the infection types have the 
same intensity, the type 2 infection can grow unboundedly from the configu- 
ration /(£), but not from where the initial disadvantage for the type 
2 infection turns out to be too severe. 

Before proceeding, we mention the following question, which has been 
pointed out to us by Itai Benjamini as well as an anonymous referee: Con- 
sider the case when the infections have the same intensity (i.e., A = 1) and 
assume that the type 1 infection occupies not only the negative xi-axis (as 
in /(£)), but a cone around it. What is the critical slope of the cone for 
which there is a positive probability for the type 2 infection at the origin 
to grow unboundedly? It seems likely, as suggested to us by Itai Benjamini, 
that the critical case is when the cone fills the whole left half-space. That the 
type 2 infection cannot survive when the type 1 infection occupies the whole 
left half-space follows from Theorem 1.1(a). That infinite type 2 growth has 
positive probability for any smaller type 1 cone remains to be proven and 
goes beyond the scope of this paper. 

The rest of the paper is organized as follows. In Section 2, some results 
are described for a one-type process started with the entire hyperplane TL 
infected at time 0. These results are then used in Sections 3, 4 and 5 to 
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prove the claims in Theorem 1.1 for A > 1, A = 1 and A < 1 respectively; 
these cases will be referred to as supercritical, critical and subcritical. Finally, 
some concluding remarks appear in Section 6. 

2. The one-type process started from Ti. In this section, we recall from 
the literature the results needed to deduce that the asymptotic speed along 
the xi-axis of the growth in a one-type process started with all sites in 7i 
infected at time is the same as the asymptotic speed along the axes in a 
one-type process started from the origin. We also state a result on the rate 
of convergence of the speed, originating with Kesten [11], and introduce a 
hampered version of the process where only edges in a certain part of Z'^ 
can be used. 

First, consider a one-type process with intensity A started from the origin. 
In such a process, each edge e of the Z'^-lattice has an exponential random 
variable T(e) with parameter A associated with it. The travel time of a path 
r is defined as 

T(r) :=^r(e) 

and, for x G Z*^, the time when the point x is infected is given by 

T^{x) := inf{r(r) : F is a path from to x}. 

Write n= (?7,, 0, . . . , 0). A basic result then, obtained from the subadditive 
ergodic theorem, is the existence of a constant fj,x> such that T^{n)/n — > 
fix almost surely and in Li as n ^ oo; see, for example, Kingman [12]. 
Defining fi = ni, it follows from a simple time-scaling argument that iJ,\ = 
fj.X~^, hence we have 

(2) lim = uA~^ a.s. and in Li. 

More generally, for any x G Z*^ with x ^0, there is a constant fJ,{x) > such 
that T^{nx)/n^ fi{x)\~^ as n — > oo. Hence, the infection grows linearly in 
each fixed direction and the asymptotic speed of growth is an almost sure 
constant. That this is also true when all directions are considered simul- 
taneously is stated in the shape theorem, which is the main result for the 
model. To formulate it, write SP{t) for the set of infected points at time t in 
a process started from the origin and let S,^{t) C M*^ be a continuum version 
of £,^{t) obtained by replacing each x G £P{t) by a unit cube centered at x. 

Theorem 2.1. There exists a compact convex set A such that, for any 
e > 0, almost surely 

Pit') 

(l-e)AAc^y^c(l+e)A^ 

for large t. 
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An "in probability" version of this result was established by Richardson 
[13] and strengthened to the present form by Kesten [13]. Note that A is the 
unit ball in the norm defined by that is, A = {x: fi{x) < 1}. 

Now, consider the Richardson model started with all sites in the hyper- 
plane Ti. infected. We will show that the asymptotic speed of growth in the 
direction of the rri-axis in such a process is, in fact, the same as the speed 
of growth along the axes in a process started with only the origin infected. 
To this end, write T^(x) for the time when the point x is infected in 
a process started from 7i. 

Proposition 2.2. In the unit rate one-type process, we have, as n —> oo, 
that r^(n)/n ^ fi in Li. 

Proof. Let Tin be the hyperplane at xi-coordinate n, that is. Tin = 
{x G Z'^ixi = n}. (In this notation, we have Ti = Tio.) The first important 
observation is that, in the first-passage percolation representation of the 
model that is used in this paper, the time when the point n is infected in 
a process started from Ti is the same as the time when the first (in time) 
point belonging to Ti is infected in a process started from n. Furthermore, the 
minimal travel time from n to has the same distribution as the minimal 
travel time from the origin to Tin, which we denote by T^{Tin)- Hence, we 
have that 



r^(n) , T^iTin) 



where = denotes equality in distribution. Just as in establishing (2), the sub- 
additive ergodic theorem can be applied to show that T^{Tin)/n converges 
in Li (and almost surely) to a constant c. Using (3), the proposition follows 
if we can show that c = ji. Clearly, since n S Tin, we have T^iTin) < 
which gives c < lim„T^(n)/n = /i. The reverse inequality follows from the 
fact that the asymptotic shape of a process started from the origin is con- 
vex. Indeed, having c < /x would contradict the convexity of the asymptotic 
shape A stipulated in Theorem 2.1. Hence, c = fi, as desired. □ 

Next, we state a result on the convergence rate in (2). The result origi- 
nates with Kesten [11] and is formulated here for the time T^{Tin) when the 
hyperplane Tin is reached by the infection. In the original formulation, the 
estimate concerns the time T^{n) when the single point n is infected, but it 
is pointed out that the bound also applies to passage times to hyperplanes. 
Kesten [11] also contains results on the convergence rate in the shape theo- 
rem; related results, improving some of the bounds of Kesten, can be found 
in Alexander [1]. 
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Theorem 2.3 (Kesten [11])- There exist constants ci, C2 and C3 such 
that 



r°(w„)-E[T«(w„)] 



>x]< c\e ^'^^ for X < c^n. 



n 

Combining this estimate with Proposition 2.2 gives the following lemma, 
which will be useful in controlling the type 1 infection in a two-type process 
started according to I{7i). 

Lemma 2.4. In a unit rate one-type process, for any e > 0, there exist 
constants c and c' such that 

(4) P(r^(n) < (1 - e)nfi) < ce'"'^ 

for large n. 

Proof. Trivially, we have 

P(r^(n) < (1 - e)n^) < P(|r^(n) - n/i| > en^i). 

To obtain |T^(n) — n/i| > en^, by the triangle inequality, at least one of the 
quantities |r''^(n) - E[r^(n)]| and |E[r^(n)] - n/i| must exceed en/i/2. By 
Proposition 2.2, we will, indeed, not have |E[r^(n)] — n//| > en/u/2 when n 
is large, hence 

P(|T^(n) - n^l > en^i) < p(^|T^(n) - E[r^(n)]| > 

for large n. As observed in the proof of Proposition 2.2, the passage time 
r'^(n) has the same distribution as T^{7in)- Choosing x = ^/nefi/2 (which 
is clearly smaller than c^n when n is large) in Theorem 2.3, we hence obtain 

P(^|r^(n) - E[r^(n)]| > ^) < 

and the lemma follows. □ 



Using Lemma 2.4, it turns out that we can establish that the convergence 
in Proposition 2.2 also holds in the almost sure sense. 

Proposition 2.5. Asn ^ 00, we have that T^(n)/n — > /i almost surely. 

Proof. We need to show that for any e > 0, 

(5) limsup — <{l+e)^i 
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almost surely and that 

T^(n) 

(6) liminf —>il-e)fi 

n — ^oo 11 

almost surely. First, note that (2) implies that (5) holds with T^{n) in place 
of T^(n). But, obviously, T^(n) <T°(n), so (5) is established. 

Next, in order to show (6), fix < oo in such a way that (4) holds for all 
n> N. From Lemma 2.4, we then have that the expected number of 77, > 
for which the event {T^(n) < (1 — e)nfi} happens is at most 

00 

n=N 

which is finite. Hence, by the Borel-Cantelli lemma, (6) is established. □ 

The last consideration before moving on to the two-type process is to 
show that the growth of a one-type process restricted to a "tube" around 
the xi-axis behaves approximately as an unrestricted process when the tube 
is large. This will be needed to control the type 2 infection at the origin 
in a two-type process started from the configuration To formulate 

the result, we introduce a new, hampered version of the one-type process 
by placing "walls" in restricting the growth in all directions except one. 
More precisely, we consider a process with the same dynamics as the original 
one, but where only sites in the set 

(7) Jli,:={xGZ'^:|3;,| <6for ahi/l} 

are susceptible to the infection. We write [t] for the set of infected points 
at time t in such a process, started with a single infection at the origin 0. 
The following lemma says that is filled with infection linearly in time. 

Lemma 2.6. Consider a hampered one-type process with rate A. For any 
dimension d, there is a real number fix^fy > such that, for any e G 
almost surely 

{x G : |xi| < (1 - e)tfi^l} C C {x G : |xi| < (1 + e)i/i^]} 

for all sufficiently large t. 

Being a completely standard adaptation of the proof of the shape theorem, 
the proof of this lemma is omitted. 

The constant ^\^b is the analog of /iA~^ in the unhampered process. That 
is, if r^*(x) denotes the time when the point x is infected in a hampered 
process, then we have fix^b = lim„T''*(n)/n. When b is large, it is reasonable 
to expect that the speed of growth in the hampered process is close to the 
speed of an unhampered process and hence that /iA,b is close to /iA~^ for 
large b. This intuition is confirmed by the next lemma. 
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Lemma 2.7. As b — > oo, we have that — > /iA . 

The proof of this lemma is again a straightforward adaptation, this time 
of the proof of Lemma 4.4 in Deijfen, Haggstrom and Bagley [4], where the 
same result is established for a continuum counterpart of the Richardson 
model. It is therefore omitted. 

3. The supercritical case. Our main task in this section is to prove co- 
existence when A > 1. 

Proposition 3.1. For any A > 1 and any d>2, we have 

Here and later, we will make use of the following convenient way of con- 
structing the two-type Richardson model: attach to each edge e in the edge 
set E^d of the Z'^-lattice independently two independent exponential vari- 
ables Ti(e) and T2(e) with mean 1 (resp. A~^), indicating the time it takes 
for the type 1 (resp. type 2) infection to traverse the edge. 

Proof of Proposition 3.1. Fix A > 1. For a hampered one-type pro- 
cess in a region Qh, as in (7), Lemma 2.7 tells us that liuih^oo ^J'X,b nX^^ 
and, since /iA~^ < /x, we have fix^i) < /x when b is sufficiently large. Fix such 
a b and define 

= {x G > 1}. 

For X G Z*^, define S'^'^\'f°J'(x) as infr J^eer ''"i(^)' where the infimum is over 
all paths starting at 7{\ {0} and ending at x. Furthermore, for x G il^, 
define S'^'^{x) as infr X^eer ''"2(e), where, this time, the infimum is over all 
paths starting at and passing through vertices in ilt, only. 

The point of these definitions is the following observation, which is easy 
to see and which turns out to be instrumental in proving Proposition 3.1. 

Suppose that at least one vertex in is eventually infected by 
type 1 and let x be the first (in time) vertex in Jl^ for which this 
happens. We then have 

Si'^\W(x)<5,''°(x). 

So, if we can show that 

(8) P(5i'^\^°>(x) > Sl'^ix) for ah xen+)> 0, 

then we know that, with positive probability, all x G 0^ eventually become 
infected by type 2, and the proposition follows. 
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To do this, define fi* = ^"'"2^''' so that < A^* < and note that, by 
Proposition 2.5, we have the existence of some (random) Afi < 00 such that 

S^''^^\W{x) > ^*xi for all x G 17+ with xi > Mi. 

Likewise, Lemma 2.6 guarantees the existence of some (again random) M2 < 
00 such that 

S'^'°(x) < /i*xi for all x G with xi > M2. 

Taking M = max{Mi, M2} yields 

S^''^^\W{x) > Sl'°{x) for all x e with xi > M, 

which is very close to proving (8). To rigorously get from here to (8), we 
employ the following conditioning argument. Choose an m such that 

(9) P(S^'^\^°>(x) > Sl'°{x) for aU x G f7+ with xi > m) > 
and define E'^'"^' as the set of edges {x,y) G E^d such that either 

x = and y = 1 = (1, 0, . . . , 0) 

or 

X, y G $7^, xi < m, yi < m + 1. 

We write Dm for the event in (9). We will condition on the ri(e) variables 
for all edges e G E^^d, together with the r2(e) variables for all edges e G 
Ej^d \ E^'"^. Let ^ be the event that 

(10) P(^m|{Ti(e)}e6i=;^,,{r2(e)}^g^^^^^+.™) >0 

and note that since P(Z)m) > 0, we must also have P(A) > 0. On the event A, 
there exists a random 7 > (depending on {Ti(e)}ee£;^^ and {^2(6)}^^^ d\E^'" 

such that if T2(e) < 7 for all e G E^'^ , then occurs. By further decreas- 
ing 7 > 0, we can also ensure (due to the fact that S^"^"^^^^ ix) > for all x 
almost surely) that if T2(e) < 7 for all e G E^'^ , then we also have 

5i,W\{0}(^^ > 5fe^'°(x) for ah x G S^^^ with xi < m. 

Finally, note that on the event A — that is, under the conditioning in (10) — 
the probability that T2(e) < 7 for all e G E^''^ is strictly positive; this is 
simply because the edge set E'^'"^ is finite and 7 > 0. Hence, 

p(5iM{0}(^) > Sl'^ix) for all x G J7+|A) > 0. 

Since P(^) > 0, (8) follows and the proof is complete. □ 

The other result concerning the supercritical case that we need to prove 
is the following. 
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Corollary 3.2. For any A > 1 and any d>2, we have 

^£,0(^2) >o. 

One way to prove this, is to note that the proof of Proposition 3.1 can 
be easily adapted to handle the corollary. Or, to be a bit more careful, we 
can invoke the proposition itself, together with the following easy lemma. 
For disjoint subsets ^1 and ^2 of Z'^, write P^^'^^ for the law of the two- type 
process with parameters Ai and A2, with all sites in ^1 initially infected by 
type 1, all in ^2 by type 2 and all others uninfected. 

Lemma 3.3. Suppose that ^1 and ^2 disjoint subsets of 7,'^, and ^'i 
and ^2 likewise. If (,1 ^ (,[ and ^2 5 ^2^ th^^ 

Proof. Couple the two processes using the same Ti(e) and T2(e) vari- 
ables. Writing (Ci(^),C2(i)) for the state of the first process at time t in 
the obvious way and similarly for the second process, it is straightforward 
to show that the relations ^i{t) C ^[{t) and ^2(*) 2 ^2(*) preserved as t 
increases. Letting t ^ 00 proves the lemma. □ 

Proof of Corollary 3.2. Define C in the same way as C, except 
with the roles of the xi- and X2-coordinates interchanged, that is, let 

C' = {x : X2 < and Xj = for all i / 2}. 

Then, by symmetry, we have Pj-^,^Q{G2) = P£'q(G2). But, since C C H, we 

can invoke Lemma 3.3 in order to deduce that P^]^q{G2) > P-ll'''o{G2) , the 
latter probability being positive by Proposition 3.1. □ 

4. The critical case. The critical case A = 1 is, in some ways, the most 
interesting, particularly in view of the fact that whether or not a single 
type-2 infection at the origin has a chance to survive against an infinite 
army of type-1 enemies depends on whether these initially get to occupy all 
oi 7i\ {0} or just £ \ {0}. The two results we need to prove in this section 
are the following. 

Proposition 4.1. For any d>2, we have -P^'o(G2) = 0. 
Proposition 4.2. For any d>2, we have P^'o(G2)>0. 
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For A = 1, it is often convenient to replace the construction following 
Proposition 3.1 by one that involves not two but just one travel time variable 
r(e) associated with each e S E^d. Here, the r(e)'s are taken to be i.i.d. 
exponentials with mean 1 and represent the time it takes for either of the 
two infection types to traverse the edge. We will employ this construction in 
the proofs of both Propositions 4.1 and 4.2. A nice feature of the construction 
is the following. Suppose that the process starts at time with the nodes in 
^1 C Z'^ infected by type 1 and the nodes in ^2 C Z'^ infected by type 2. We 
write T^^{x) for infr X)eer ''"(^)> where the infimum is over all paths starting 
in ^1 and ending at x, and define T^^{x) analogously. Then, x becomes 
infected precisely at time mm{T^^ {x),T^^ (x)} and, furthermore, 

(11) X gets infected by type 1 if and only if T^^(x) < T^^{x). 

The following lemma (an easy variation of Lemma 3.3) will be useful in 
the proof of Proposition 4.1. 

Lemma 4.3. Consider the (one-type) Richardson model with parameter 
1 on a bounded degree graph L with vertex set V and edge set E starting at 
time with the set ^ C ^ initially infected. Also, consider the same model 
with parameter 1 on another graph L' = {V',E') starting at time with 
the set S^' C V initially infected. If L is a subgraph of L' , in the sense that 
y C y and E <^ E' , and, furthermore, we have ^ C then, for any t > 
and any t] C E, the probability that all vertices in rj are infected at time t is 
no greater for the process on L than for the process on L' . 

Proof. Simply couple the two processes in such a way that for each 
eE E, they use the same T(e) variable. □ 

Proof of Proposition 4.1. We wiU show that 
P(type 2 infects infinitely many sites in the half-space {x : xi > 1}) = 0. 

(12) 

Once that is done, we can, by symmetry, infer the corresponding statement 
for the other half-space and thus conclude that P(G2) = 0. 

Recall, for integer 6, that Tib is the set of vertices x G Z'^ whose xi- 
coordinate is b (so that, in particular, TIq = Ti). For x G and y Ti, 
we write y — > x for the event that the infimum T^(x) is obtained by a path 
from y to X. Beginning with all nodes in TC infected, we think of y ^ x as 
meaning that the infection eventually hitting x descends from y. Since the 
r(e)'s are independent with a continuous distribution, we have, for fixed x, 
that the event y ^ x happens for only one y almost surely (i.e., there are 
no ties). Thus, 

(13) 5]P(y^x) = l. 

yen 
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Fix 6 > 1, y,y' £7i and x, x' € Tib in such a way that x and y differ only 
in their first coordinate, likewise for x' and y' . Symmetry implies that 

P(y ^ x') = P{y' ^ x) 

and, in conjunction with (13), this implies that 

(14) ^P(2/^x) = l 

xeHb 

for any y EH. Let Xf, denote the number of sites x £Ti.b such that — > x. 
Then, to prove (12) is the same as showing that 

(15) p|^^Xfe = cx)^ =0. 

By (14), we have E[X6] = 1, whence P{Xf, = oo) = for any b, so in order 
to prove (15), it is enough to show that 

(16) f( Urn Xb = o] =1. 

Ideally, since E[Xb] = 1 for each b, we would now like to endow the sequence 
{Xi,X2, ■ ■ ■) with a martingale structure (with respect to some filtration). 
Since a nonnegative martingale converges almost surely and this one pre- 
sumably could not converge to anything other than 0, that would settle 

(16) . But we cannot see how to do this (it is probably not even possible) 
and will settle for a different solution that, although a bit less clean, is still 
reminiscent of a martingale approach. 

To this end, we first define, for x G with xi <b and y ^Tl, the passage 
time T^{x) as infr X^eer where the infimum is over all paths from y to 
X that do not pass through any vertex z with z\ > h. Also, define 

(17) r,«(x)=infT^(x) 

and for y' G let y' ^ x denote the event that the infimum in (17) is 
attained for y = y' . Finally, write X^ for the number of vertices in fib such 

that x. By the same argument as for Xb, we have, for any 5 > 1, that 

(18) B[X*b] = 1. 
We now claim, crucially, that 

(19) {Xb > 1} implies {X^ > 1}. 

To see this, assume that Xb > 1 and take x to be the vertex among those 
in Tib satisfying ^ x for which T^{x) is smallest. Then, the path from 
to X cannot pass through any other vertex z in Tib (because if that were 



THE TWO-TYPE RICHARDSON MODEL 



13 



the case, z would satisfy — > z with a smaller T°). Hence, T^{x) = T°(x). 

Since Tfj{x) > Ty{x) for all y £Ti., it follows that x and the claim (19) 
is warranted. 

Now, write G2 for the event that limsupj,^oo > 0. Using (19), we can 
show that (16) — and thereby the proposition — follows if we can show that 

(20) P(G'^) = 0, 

so this is what we set out to prove. 

Our argument will involve the filtration {^bj^Q, where J-'b is the u-field 
generated by the r(e)-variables for all edges e = {x,y) with xi,yi < b. Note 

that, for X G Tib, tbe event x and the random variables T^(x), T^{x) 
and are all .T^j-measurable. Levy's 0-1 law (see, e.g., Durrett [5]) tells us 
that 

(21) limP(G^|.;^5) = lG* 

almost surely, where Ig* is the indicator of the event . The conclusion of 
the argument will be to show that, with probability 1, 'P{G2\^h) does not 
converge to 1. 

Fix e > small, in such a way that n='^ is an integer. For x,?/ S Z*^, 
write dist(x,y) for the graph-theoretic distance between x and y in the Z'^ 
lattice. We claim that there exists a n < 00 such that, for any 6, the event 
Db,u has probability at most |, where we define Db^u ^ the event that 

< n, while, for some x,y ^TLb such that x and dist(x, y) < 
n, we have Tl^{y)>Tl^{x)+u. 

To see this, we will invoke a comparison with the classical one-type 
Richardson model on Z*^"^ with infection rate 1 and starting at time with 
a single infection at the origin. Since, almost surely, any finite set of sites is 
eventually infected in this model, we can find a u such that the probability 
that all sites within distance n from the origin are infected by time u is at 
least 1-^. 

Returning to the model in Z'^, consider (for some fixed A; G {1, . . . , n}) the 
vertex x that has the kth smallest value of T^{x) among those for which 

X (provided at least k such vertices exist). For a fixed realization of 
the process up to time T^{x), a stochastic domination argument now shows 
that the probability that all sites in Tib within distance n from x are infected 
by time T^{x) -|- u is at least 1 — to see this, apply Lemma 4.3 with 
L equal to Tib with edges between Euclidean nearest neighbors (this graph 
is isomorphic to Z*^"^) and ^ = {x}, L' equal to Z'^ restricted to vertices z 
with zi <b and ^' equal to the set of vertices infected by time (x) in the 
original process. 
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Summing the resulting complementary probability bound ^ from 1 to n 
gives the desired bound 

P(A,J<|. 

Furthermore, Markov's inequality applied to (18) yields P(^^ > n) < ^ = ^. 
If we now define the event 

Di,,u = Db,u^{X;>n}, 

we obtain 

P(56,„)<| + | = e. 

Let -1 denote set complement. The next crucial claim is that (given e > 
as above) there exists a 5 > (independent of b) such that the event -'-Dfe,^ 
implies 

(22) p{x;^, = o\j^k)>s. 

To see this, assume that ^Di, y_ occurs and consider the following event Ai,^u 
which, together with -iD^ is enough to guarantee that X^^^ = 0. Namely, 
let Ab^u be the event that: 

b 

• for all e = {x,y) with y G TLb+i, x G Tib, -/^ x and x within distance n 

from some z G Tih such that z, we have T(e) < 1; 

• for all e = (x, y) with x, ?/ G 'Hh+i within distance n + 1 from some z G Tii, 

such that z, we have r(e) < 1; 

• for all e = with y G TLb+i, x £7ib and x, we have r(e) > u + n+1. 

The point of this definition is that if ^D^^u and Ai, y_ occur, then none 

of the vertices z G Tib with z will have time to infect their neighbor in 
"Hb+i before infection creeps in from another direction, and therefore X^^-^^ 
will equal 0. These requirements on Ab^u are guide restrictive, but note that, 
on the event -^Db^u^ there are at most n edges e for which we require that 
r(e) > ti + ?i + 1 and there are at most nd{2n)'^~^ edges e for which we 
require that r(e) < 1. Thus, on the event -^Db^u, we have 

P{Ab,M > (e-(-+"+i))"(l - e-i)"'^(2n)''-\ 

which is a small number indeed, but strictly positive, and the whole point of 
this exercise is that it does not depend on b. Thus, we have for any b on the 
event ^Db,u that (22) holds with 5 = (e-("+"+i))"(l - e-i)""'(2n)<^-i ^ g^^^^g 
X^_^^^ = precludes G2, we also have on the event -^Db^u that 

P(G^|^6)<l-5. 
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So, in order for the limit limf,^oo P(G2|^;,) in (21) to equal 1, we need Db^^ 
to occur for all sufficiently large b. But Db^^ was defined in such a way as to 
guarantee that F{Db^u) ^ ^ for any b, so we can conclude that P(G2) < £■ 
Since e > could be taken to be arbitrarily small, we obtain P(G2) = and 
the proof is complete. □ 

That was rather involved. Fortunately, the proof of Proposition 4.2 is 
somewhat more straightforward. 

Proof of Proposition 4.2. Recalling the notation n= (n,0, ...,0), 
the initial configuration in Proposition 4.2 consists of a single type-2 in- 
fection at 0, competing against type-1 infections at —1,-2, Intuitively, 

the best hope for type 2 seems to be to rush off along the positive xi-axis. 
Assuming this to be the case and viewing the model in terms of the T(e)- 
variables as before, we will set out to prove that 

(23) liminfpf T^(n) = inf r'"(n)^ > 0. 

n^oD \^ m<0 / 

Note that the event in (23) is precisely the event that the vertex n becomes 
infected by type 2. Hence, if we can prove (23), we can deduce that, with pos- 
itive probability, infinitely many vertices on the positive xi-axis are infected 
and the proposition will follow. 

Symmetry implies that for any n > 1, we have 

(24) P f r°(n) = inf T"^(n)) = P f r°(n) = inf r°(l) 

\ m<0 / \ l>n 

We will work with the right-hand side of (24), the advantage of this being 
that it has a useful interpretation in terms of the one-type Richardson model. 
Namely, for the one-type model starting at time with a single infected site 
at the origin, the event {r''(n) = inf/>„T^(l)} is precisely the event that 
the node n is infected before any node further away on the positive xi-axis 
is infected. 

To this end, for the one-type Richardson model, define Y{t) as the number 
of nodes on the positive xi-axis that are infected by time t. Further, define 
Y^{t) as the number of nodes on the positive xi-axis that are infected by 
time t, with the additional property that at the time of their infection, they 
became the rightmost infected node on the xi-axis. 

We have, from (2), that 

t—^00 t 

almost surely, where // > is the time constant discussed in Section 2. What 

Y — ^ (I) 

about — At any time t, there exists a rightmost infected node on the xi- 
axis and this node infects its neighbor-to-the-right at rate 1. Every time such 
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an infection occurs, Y'^{t) increases by 1. Hence, the process {Y^ (t)}t>Q 
stochastically dominates a rate-1 Poisson process, whence 

(25) lim inf ^>1 

t^oo t 

almost surely. Again by (2), we have, with probability 1, that for any e > 0, 
eventually all infected nodes on the positive xi-axis at time t have an xi- 
coordinate that does not exceed (1 + e)t/x~^, that is, the number of infected 
nodes on the positive xi-axis at time t does not exceed (1 + e)tfi~^. In 
conjunction with (25), this implies that 

n 

lim^inf ^ I{TO(j)=i„f,>^. ro(i)} > 
i=i 

almost surely. Hence, by Fatou's lemma, we have 

lim^infn~^^^p(^r°(j) =infr°(l)^ > ^t, 

implying that 

lim inf pf r°(n) = inf r°(l)^ > u. 

n^oo \ l>n ^ ' ) 

Using the identity (24), this implies (23), and the proof is complete. □ 

5. The subcritical case. For the subcritical case A < 1, we shall see that 
the following result holds, which is more general than the subcritical cases 
of both part (a) and part (b) of Theorem 1.1. 

Proposition 5.1. Consider the two-type Richardson model on , d> 
2, with types 1 and 2 having respective intensities 1 and X, starting with 
type 1 in d and type 2 in ,^2 C Z*^. // ^1 is infinite, ^2 is finite and 
A < 1, then the event G2 of unbounded survival for type 2 has probability 0. 

It turns out that this result is a direct consequence of the following propo- 
sition which was instrumental in proving the main result in Haggstrom and 
Pemantle [8] mentioned in Section 1. Let A be the asymptotic shape for the 
Richardson model, as defined by Theorem 2.1. 

Proposition 5.2. For any A < 1 and any e > 0, we have 

lim supP3;'^. (G2) = 0, 

where the supremum is over all initial configurations (^1,^2) such that 
^2 is contained in rA, while 

(26) 

^1 is not contained in (1 -\- e)rA. 
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In fact, the proposition as stated in Haggstrom and Pemantle [8] dealt 
only with the case where was finite, but the generalization to infinite 
^1 follows immediately from Lemma 3.3. Proposition 5.1 now follows from 
Proposition 5.2 upon noting that if is infinite and ^2 is finite, then the 
pair (^1,^2) satisfies (26) for all sufficiently large r. 

6. Concluding remarks. To see that our main result. Theorem 1.1, has 
now been proven is just a matter of collecting the results from the previous 
three sections. Theorem 1.1(a) follows from Propositions 3.1 (supercritical 
case), 4.1 (critical case) and 5.1 (subcritical case), while Theorem 1.1(b) 
follows from Corollary 3.2 (supercritical case) and Propositions 4.2 (critical 
case) and 5.1 (subcritical case). 

We end the paper with the observation that Proposition 4.2 allows us to 
construct a simple proof of the fact that infinite coexistence is possible in the 
critical case A = 1 starting from finitely many infected nodes of each type 
(recall the result of Deijfen and Haggstrom [3] that the particular choice 
of finite initial configurations does not matter, as long as type 1 does not 
already "strangle" type 2 or vice versa). There already exist several proofs 
of this result — Haggstrom and Pemantle [7] produced one for d = 2, while 
Caret and Marchand [6] and Hoffman [9] proved the result for arbitrary 
d — but, since the result is central to the study of the two- type Richardson 
model, we feel it is worth the extra effort to state a new, simple proof. 

Theorem 6.1. For the critical (X=l) two-type Richardson model on 
'L^ in any dimension d>2, there exists an n such that if the model starts 
with a single type-1 infection at and a single type-2 infection at n, then 
infinite coexistence has positive probability. 

Proof. Consider the usual edge representation of the critical model, 
where each edge is assigned a r(e) representing the time it takes either 
infection to traverse it. Also, as usual, for ^ C Z*^ and x G Z*^, write T^{x) 
for the sum of the r(e)'s along the fastest path starting in ^ and ending at 
X. Define two random sequences {Xn}^=_oo and {^n}^_oo ^ follows. Set 



and 



Yn 



1, if r"(z) < T^-,'^-^,n-2,n-i} 1^^^ infinitely many z G Z^, 
0, otherwise, 



1, if r"(z) < Ti"+i'"+2.»+3v}(2) for infinitely many z € : 
0, otherwise. 



Let a = P{Xq = 1). Proposition 4.2 tells us that a > 0. The process {Xn} 
is stationary, so F{Xn = 1) = a for any n. By symmetry, P(Yn = 1) = a also 
holds for any n. Furthermore, {Xn}'^^_^ arises in a stationary way from 
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an i.i.d. process and is therefore ergodic, so P(X„ = 1 for some n > 1) = 1. 
Hence, we can find an n such that 



On the event {Yq = = 1}, we have (by definition of the two processes) 
that r°(z) < T"(z) for infinitely many z G Z'^ and that r"(z) < T°{z) for 
infinitely many z G Z'^. Thus, (27) guarantees that infinite coexistence has 
positive probability for the two-type model starting at and at n. □ 
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